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We investigate the quantum optical properties of a single photon emitter coupled to a finite-size 
metal nanoparticle using a photon Green function technique that rigorously quantizes the electro- 
magnetic fields. We first obtain pronounced Purcell factors and photonic Lamb shifts for both a 
7-nm and 20-nm radius metal nanoparticle, without adopting a dipole approximation. We then 
consider a quantum-dot photon emitter positioned sufficiently near to the metal nanoparticle so 
that the strong coupling regime is possible. Accounting for non-dipole interactions, quenching, and 
photon transport from the dot to the detector, we demonstrate that the strong coupling regime 
should be observable in the far-field spontaneous emission spectrum, even at room temperature. 
The emission spectra show that the usual vacuum Rabi doublet becomes a rich spectral triplet or 
quartet with two of the four peaks anticrossing, which survives in spite of significant non-radiative 
decays. We discuss the emitted light spectrum and the effects of quenching for two different dipole 
polarizations. 



PACS numbers: 42.50.Pq, 78.67.Bf, 73.20.Mf 

I. INTRODUCTION 

The route to photonic vacuum engineering tradition- 
ally employs a lossless dielectric cavity system, exploiting 
an optical mode with a suitably large quality factor, Q, 
and small effective mode volume, V . The local photon 
density of states (LDOS) scales proportionally with the 
Q/V factor. Enhancing the LDOS through the use of 
small cavities^ has shown to be a very effective method 
for increasing the radiative decay rate of an emitter via 
the Purcell effect 2 -. In solid state materials, cavities are 
created using various structural designs, including pho- 
tonic crystal lattices with defects^, and etched micro- 
pillars made of Bragg reflectors^. These dielectric cav- 
ities have shown some remarkable successes in quantum 
optics, but the lower limit on V in such systems is typi- 
cally set by diffraction, with V « (A/n) 3 , where n is the 
refractive index of the cavity. Additionally, when one 
uses quantum dots (QDs), the narrowband resonance as- 
sociated with high Q requires very long non-radiative ex- 
citon decay times, only achievable at low temperatures. 

In an effort to further increase the LDOS and decrease 
the system size to sub- wavelength dimensions, it can be 
advantageous to examine plasmonic systems where light 
is confined to the surface of a metal and decays evanes- 
cently from its surface. For example, a metal nanoparti- 
cle (MNP) supports localized surface plasmons (LSPs)£ 
that are tightly confined spatially and not limited by 
diffraction. This allows coupling between single photon 
emitters and MNPs£ which can enhance the LDOS in a 
system as small as 10-20 nm 3 . When the LDOS becomes 
large enough, it may also be possible to enter the strong 
coupling regime where instead of the irreversible process 
of the emitter decaying and emitting a photon into the 
environment (i.e., weak coupling), the emitter can re- 
versibly exchange the photon with the environment — a 



process known as vacuum Rabi oscillations. In order for 
this to happen, the coupling between the emitter and 
the environment must exceed all possible decay channels. 
Classical predictions of strong coupling behaviour have 
been made in the context of metallic dimersi, though it is 
not known if the splitting survives in the observable spon- 
taneous emission spectrum. This reversible exchange of 
energy is fundamentally interesting and can possibly be 
useful for applications in coherent control^, quantum in- 
formation processing^, and lasing/spasingi^— . With 
regards to a quantum theory of the light-matter pro- 
cesses in the strong coupling regime, several compli- 
cations arise in the theoretical description of coupling 
quantized light to a MNP, including the need to quan- 
tize the fields in a dissipative/lossy medium. Waks and 
Sridharani^ introduced a useful quantized cavity-QED 
(quantum electrodynamics) treatment of a coupled MNP 
and a dipole emitter [e.g., a QD], but the MNP was 
treated within the dipole approximation^; however, it 
is now well known that the dipolc-approximation can 
yield poor agreement with exact (i.e., non-dipole) cal- 
culations obtained from finite-size MNPs — unless placed 
a few radii from the MNP surface^—. Triigler and Ho- 
henester— have examined the strong coupling dynamics 
between a molecule and a cigar-shaped MNP employ- 
ing a mode expansion technique which incorporates the 
higher order plasmon modes; their quantum approach 
assumes a Lorentzian form for the broadening of the 
modes, via Lindblad superoperators in a master equation 
formalism~*2S; this useful non-dipole study predicts the 
strong coupling regime is possible between a MNP and 
a molecule though there is no connection to the emis- 
sion spectrum. For dielectric cavity systems, the effects 
of propagation to a detector is generally assumed to not 
change the spectral shape of the emitted photons. How- 
ever for metallic system, because of the losses associated 



with the MNP and quenching, it is important to compute 
the generalized light spectrum (i.e., away from the QD 
position) to first realize if the strong coupling effects are 
observable, and secondly, to see how the spectral signa- 
tures may change and how they would be measured. 

In this work, we develop a theoretical formalism that 
allows one to obtain the emission spectra at any spa- 
tial position of the detector. In Sec. HH we describe an 
exact medium-independent quantum optics approach — 
formulated in terms of photonic Green functions — to de- 
scribe the cavity-QED interactions and photon transport 
between a dipole emitter (QD), a finite-size MNP and a 
detector. A schematic of the nanoscale interaction ge- 
ometry is shown in Fig. (Ua). I n Sec. EE we present 
various numerical results and calculations for the cou- 
pled QD-MNP system. We first calculate the classical 
Green function above a MNP using a well established 



scattering approac 



21,22 



and subsequently calculate the 



LDOS and photonic Lamb shift from a nearby dipole 
emitter, using two different size MNPs (7 nm and 20 nm 
radius). We find significant enhancements in the LDOS 
near the MNP surface* 6 and simultaneously observe enor- 
mous Lamb shifts. We then examine the spectral proper- 
ties of a QD dipole emitter in the strong coupling regime. 
We compute the far-field spontaneous emission spectrum, 
fully accounting for non-Markovian relaxation and prop- 
agation effects to the detector. The spontaneous emis- 
sion spectrum is shown to yield clear signatures of the 
strong coupling regime, but is found to be much richer 
than the usual vacuum Rabi splitting known from simpler 
cavity-QED systems (e.g., using dielectric cavities) due 
to the interplay between higher-order mode coupling and 
dipolar-mode coupling; the non-Markovian spectra yield 
a spectral triplet or even a quartet of resonances, where 
two of the peaks anti-cross, thus signalling the strong 
coupling regime. We present the strong coupling spectra 
for two different QD-dipole polarizations and discuss the 
effects of quenching. In Sec. lIVl we give a brief discussion 
about possible experimental configurations for observing 
our predictions, and in Sec. [V] we conclude. 



II. THEORY 

A. Green function of a spherical metal-nanoparticle 

The classical photon Green function in a medium with 
e(r,u>) (complex dielectric constant) and /i = I, is de- 
scribed through the following equation: 



V x V x G(r,r'; w) - e(r, cJ)fc 2 G(r, r'; w) = k%5(r-r') , 



(1) 



where fco = u) /c, where w is the angular frequency and c is 
the speed of light. The dipole-response function (Green 
function), G, can connect to to both classical and quan- 
tum light-matter interactions. For the MNP problem of 
interest, we will discuss the Green function within and 




FIG. 1. (color online) (a) Schematic of the MNP embed- 
ded in a background material with permittivity of £(,. The 
MNP, with radius a, and permittivity, e m , is located at the 
origin. The single photon emitter (quantum dot) at r^ is lo- 
cated at height, h, above the surface of the MNP. We also 
consider a point-like detector at r_o located along the same 
axis at height ho above the metal surface, (b) LDOS peak 
as a function of height above a 20-nm (blue-dark) and 7-nm 
(orange-light) spherical silver MNP for an a;-oriented dipole. 
The non-dipole result (for the MNP) is given by the solid line 
and the dipole-approximation result is given by the dashed 
line. For comparison, selected finite-difference time-domain 
results are shown as squares for 20 nm MNPs. 



outside the dipole approximation. Typically for small 
MNPs (uiyfsba/c <C 1) of permittivity e m embedded in a 
material with permittivity e b , the Green function is ob- 
tained through the Dyson equation where we assume that 
the spherical MNP response can be modelled through the 
metal polarizability function: 



M = 



67rc ;1 a 3 J 



(2) 



with the bare polarizability (i.e., without photon cou- 
pling to the environment), 



a° m (uj) = 4ire b a 



3 [e m (uj) - e b ) 



(3) 



(e m (uj) + 2e b ) ' 

which also accounts for radiation reaction^. Considering 
the MNP to be located at position r m , then the MNP- 
dipole Green function is obtained throug h 24 ' 25 : 

G(r, r') = G (r, r') + G (r, r m ) • a m G (r m , r') . (4) 

To account for the finite-size nature of the MNP, we 
also compute the Green function outside the dipole ap- 
proximation. For these calculations we use an estab- 
lished analytical approach where the Green function is 
expanded in spherical vector functions and the boundary 
conditions are satisfied at the edge of the spher o 21 ' 22 ; we 
relegate the details of this approach to the Appendix. 
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B. Classical light-matter interactions 

An integral solution for the classical electric field can 
be written as 



E (r, w) = E° (r, w) + J dr' G (r, r'; u) ■ P (r', w) 



(5) 



where P is a polarization source. As we will show be- 
low, in quantum optics, the E and P fields become oper- 
ators, but G remains the sam o 26 ' 27 . For a dipole emitter 
at position r d , then E (r) = E° (r) + G (r, r d ) ■ a d E (r d ), 
where the dipole polarizability of the QD exciton is given 
by 



2u d d 2 /he 



{w d - lg 1 - i~/ d uj)he ' 



(6) 



with ui d the transition frequency, 7^ the non-radiative 
broadening of the QD exciton, d the optical dipole mo- 
ment, h is Planck's constant by 2ir and £0 is the permit- 
tivity of free space. Assuming a QD dipole of the form 
d = drii, then the (projected) LDOS becomes 



Im[n a ■ G (r d ,r d ;oj) ■ m] 



(7) 



where Gn = 



= KQ^eb/6n. The units of Eq. ([7]) are conve- 
niently chosen so that the LDOS is equal to the Purcell 
factor—, which describes — in a weak coupling regime — the 
spontaneous emission rate, 



7(r<*;w) 



2d 2 p(r d ;uj)G Q 

hs 



(8) 



This total electromagnetic (EM) decay rate includes both 
radiative and non-radiative coupling with the lossy envi- 
ronment; this modified decay rate obviously depends on 
the G of the medium. In order to describe photon prop- 
agation from the QD to a detector (e.g., to the far field), 
we also consider the non-local propagator, which is de- 
fined through the two space-point Green function, 



Pt] (r,r» = 



\r\i ■ G (r,r';o;) ■ n 3 ] 
Go 



(9) 



The photonic (or anomalous) Lamb shift is also con- 
nected to the Green function, and is obtained fro m 26 ' 28 



Aw (r d ;uj) 



Re[d ■ G (r d ,r d ;uj) ■ d] 



(10) 



For the Green function used in Eq. (|T0j) . i.e., with the 
same two spatial arguments Re[G(r, r)], we will neglect 
the (divergent) homogeneous-medium contribution since 
its effect can be absorbed into the definition of uj d 25 i 29 . 

The quantities introduced above (e.g., the photon de- 
cay rate and the Lamb shift) are well known, and are 
perturbative in nature (in their definition). However, 
this is not a model restriction. Indeed, the theory above 



can fully include nonperturbative light-matter interac- 
tions. To reach the strong coupling regime of cavity- 
QED, where light-matter interactions must be included 
to all orders, one requires the dipole-medium coupling 
rate, g, to be larger than any dissipation channel o 7 ' 18 . 
For a single quasimode of the metal, e.g., f m (r), g = 
g m = y /, u! m /2heo d ■ f (r<j) so that the vacuum Rabi split- 
ting, 2g « v/7em(p)7lsp/2 > 7d,7LSP- Here 7em(p) 
accounts for all EM decay processes and 7lsp is the ef- 
fective linewidth of the LSP dipole mode; for the purpose 
of the scaling argument above, we are also tacitly assum- 
ing a Lorentzian lineshape for p(u). 



C. Quantum light-matter interactions and the 
emission spectrum 

To describe the quantum light-matter interaction, we 
adopt a medium-dependent quantization procedure for 
calculating the emission spectrum from a two-level atom 
in a lossy, non-homogeneous environmen t 26 ' 2 7 . We begin 
with the complete Hamiltonian of the coupled system, 

H =fiLo d & + &~ + h J dr J duJiujiP(r,uji) ■ f(r,cj/) 
- [a+d + a-d] -E(r d ), (11) 

where (T + , a~ are the Pauli operators of the QD exciton 
(located at position r^), and f, ft are the bosonic con- 
tinuum field creation/annihilation operators of the to- 
tal electric field operator (including interactions with the 
QD), which are indexed in the Hamiltonian with contin- 
uous eigenfrequencies u)i. The electric field operator is 
related to the bosonic field operators through 28 , 



E(r, t) = E°(r, t) + J duji J dr'G(r, r': 



u>i) 



■ v /ej(r / ,wi)f(r / ,wi;t)+H.c., (12) 

where ej is the imaginary part of the permittivity and 
E°(r,t) is the free field, i.e., the field that exists with- 
out the presence of the QD. To proceed we will adopt 
the weak excitation approximation, so that we assume 
at most one quantum in the system (this approximation 
is exact when the initial field is in vacuum). Using the 
Heisenberg equations of motion, and Laplace transform- 
ing to the spectral domain, we can subsequently obtain 
explicit expressions for <t + ,<t~ and f , ft . The total elec- 
tric field operator is ther* 2 ^: 

E(r,w) = E (r,w) 

ImG (r,r d ;uji) ■ d <7~(w) + ct + (lo) 



7T£ 
1 



w - U>1 



E°(r,t) + -G (r,r d ;u>) ■ d[a~ (lj) + a+ {u)], 

(13) 
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FIG. 2. (color online) (a) LDOS versus frequency 2 nm above a 20-nm (blue-dark) and 7-nm (orange-light) spherical silver 
MNP. The non-dipole (exact) result is given by the solid line and the dipole-approximation result is given by the dashed line, 
(b) pSi(rD, r<j) versus frequency for hd — 2 nm and hn = 1 fJ-m, (c) Lamb shift versus frequency. In (a)-(c) we use an a;-oriented 
dipole. (d)-(f) as (a)-(c) but for a z-oriented dipole. For clarity, dipole results are multiplied by a factor of 50 (20 nm particle) 
and 4 (7 nm particle) in graphs (a),(c) and (d),(f). 



in which we have used the relation % -r. — = 

ttS (cji — u>) + iP ( jjp^j ) has been used, with P is the 
principle value. The dipole operators are given by, 

<r~(ui) + <T + (a;) = 

-j [&-(t = Q)(lj +uj d ) + £+(t = Q)(a; - cj d )] 
u 2 d - ui 2 - 2w d d-G(r d ,r d ;w) • d/he 

The light spectrum is defined through S^r, ia) = 
/ °°dti / ° o ^ 2 e^( t2 -* 1 )([E(r,t 1 )]tE(r,t 2 )), which gives 

S(r,w) = ([E (r,£j)]TE (r, uj)). For r = ro, and assuming 
an initially-excited QD exciton in vacuum, one obtains^ 
the emitted light-spectrum, analytically, 

S(t d ,u) = 

d ■ G (r D , r d ; us) (us + us d ) /e 2 

u> 2 d - ui 2 - iu)^ d - 2u) d d • G (r d , r d ; uj) ■ d/Tie 

(15) 

where the point detector is assumed to be at position ro 
above the center of the MNP. We highlight that this fi- 
nal spectrum is exact in both weak and strong coupling 
limits. In order to more clearly extract the physics as- 
sociated with propagation and quenching, we will also 
examine the dipole or polarization spectrum: 

P(a;) = <*+(a;)*-M) 



ojj - uj 2 ~iu)"f d - 2u d d G (r d , r d ; u>)-d/h£ 

(16) 



which contains important information about the local dot 
dynamics. Worth to note is that Eqs. (fT5|) and (IT6l) arc 
applicable in any lossy, non-magnetic inhomogeneous sys- 
tem, provided it is possible to calculate the Green func- 
tion, which illustrates the strength of our technique. We 
also remark that it is relatively straightforward to include 
multiple QDs within this formalism^,. 



Before closing this theory section, we make a few gen- 
eral comments on the form of the QD non-radiative decay 
rate, jd- This broadening mechanism is likely caused by 
electron-phonon scattering and pure dephasing, which is 
especially important at elevated temperatures. Although 
we have not distinguished the mechanism of pure dephas- 
ing from an effective decay rate in the polarizability, the 
computed spectrum maintains precisely the same spec- 
tral shape for our chosen initial conditions^; so the dis- 
tinction of pure dephasing is not necessary for computing 
the vacuum spectrum. However, if one knows the precise 
spectral form of the QD polarizability, including the in- 
fluence of electron-phonon scattering, then only a small 
modification is needed in the above formulas^. For the 
calculations that follow below, we will adopt broadening 
values similar to colloidal dots at room temperature 32 , 
with 7^ = 10-20 meV. Note also, that since the dominant 
decay is from non-radiative coupling to the lossy MNP, 
the details of the bare exciton decay are less important 
here (e.g., in comparison to coupling to a dielectric cavity 
system) . An alternative quantum optics approach can in- 
clude phonon interactions at the level of a polaron master 
equatio n 31 ' 34 . 
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FIG. 3. (color online) (a) + (c) 7-nm particle [(b) + (d) 20- 
nm particle] with an emitter 2 nm from the surface. Graphs 
(a)-(b) and (c)-(d) show the normalized effective particle 
and far-field spontaneous emission spectra, respectively, us- 
ing the non-dipole result (blue-dark line) and the dipole- 
approximation (orange-light line) for an a;-oriented dipole. 
Transition frequencies are indicated by black dots on the 
curves. The thin grey line in all graphs indicates the LSP 
resonance (at the maximum of a m ). 



III. RESULTS 



Weak coupling regime: Purcell factors and 
Lamb shifts 



For our numerical calculations, we assume a MNP 
with a permittivity given by the Drude model, e m = 
Eoo — us^/ (lu 2 — ij m uj), and take the parameters typical 
for silver—: = 6, u m = 7.90 eV and 7 m = 51 meV; this 
gives an estimated jlsp = 60 meV and ^lsp = 75 meV 
for 7 nm and 20 nm particles, respectively, in the regime 
where the dipole approximation is valid. We consider a 
dipole emitter located h = 2 nm above a 7-nm and a 20- 
nm MNP. For the single photon emitter (QD exciton), 
we consider both x-oricntcd and z-oriented dipoles with 
a dipole moment of d = 24Debye (« 0.5 enm) which is 
comparable to (or less) than the dipole moment used in 
other works that model QDs coupled to metals^^. 

In Fig. [1] (b) we show the LDOS versus height using 
both the non-dipole and dipole calculations. We observe 
convergence between the solutions with the analytic and 
the dipole-approximation only for h > 2a, in agreement 
with Ref. HE Additionally for the 20-nm radius MNP, 
we plot the same calculations performed using finite- 
difference time-domain calculations^ (squares) and a 1- 
nm grid size (finite-size emitter); we observe excellent 
agreement between these two different methods. 

In Figs. HJa) and Hd) we plot the LDOS MS iX func- 
tion of frequency for h = 2 nm above the 7-nm and 
20-nm MNPs for rr-oriented and z-oricntcd dipoles, re- 
spectively. We immediately notice that the LDOS peaks 



FIG. 4. (color online) As in Fig. (3), but with a z-oriented 
QD dipole. 



are far separated in energy when compared to the dipole 
result, which is caused by the essential contribution from 
higher-order modest. We also see that the LDOS peak 
for both nm-size particles is comparable, but the LDOS 
peak is slightly shifted between the two different sized 
particles. When comparing between x-oriented and z- 
oriented dipoles we see that the LDOS is larger for the 
latter case by about a factor of 2. Also note that the 
difference in the dipole-approximation for the 7-nm par- 
ticle compared to the 20-nm particle is mainly due to the 
fact that they have different center to center distances; 
the shorter distance gives a larger result because of the 
scaling of the free space Green function in the near field, 
i.e., Gf ree (rd, I'm) r m | ^. 

We next consider the non-local propagator in Figs. [2] 
(b) and (e) for x-oriented and z-oriented, respectively; 
this propagator is needed to account for light propaga- 
tion from the dipole emitter to the detector. The detec- 
tor is assumed to be at a height of 1 /im above the MNP 
surface. For the 20-nm particle, the non-dipole calcula- 
tions for p nl is spectrally peaked near ui w 2.76 eV, which 
does not coincide with the peak of the LDOS (« 2.97 eV); 
however, the peak in the p nl using a dipole-approximation 
is shifted to 2.79 eV. We also observe an additional peak 
located near 2.9 eV, and we show below how this complex 
lincshapc affects the spontaneous emission spectrum. We 
can contrast these 20-nm MNP findings with the 7 nm 
results, where the p nl in the dipole-approximation agrees 
quite well with the exact result — although we begin to ob- 
serve a small shoulder in the non-dipole result which in- 
dicates a second peak. Both non-local propagator peaks 
in this region are located at 2.8 eV, which is the same 
location as the dipole peaks seen in the LDOS and again 
the difference between x-oriented and z-orientcd dipoles 
is about a factor of 2; but now the xa;-component of 
the non-local propagator is the larger (suggesting less 
quenching) . 

Figures^] (c) and (f) show the photonic Lamb shift for 
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FIG. 5. (color online) (a) 7-nm particle [(d) 20-nm particle] with an emitter 2 nm from the surface tracking the particle spectral 
peaks (red circles) and light emission spectral peaks (blue crosses) as a function of QD frequency for an i-oriented QD. The 
dashed lines correspond to the location of the QD transition frequencies used in the graphs in Fig. [3] (b) Integrated particle 
spectra for an s-oriented QD normalized by the integrated particle spectrum without the MNP (i.e., free space) as a function 
of QD transition frequency for 20 nm MNP (blue-dark line) and 7 nm MNP (orange-light line) for dipole (dashed line) and 
non-dipole (solid line) calculations, (c) Integrated particle spectrum 2 nm above a metallic half space normalized by integrated 
particle spectrum in free space as a function of QD transition frequency for z-oriented QD (blue-dark line) and i-oriented QD 
(orange- light line), (e) as (b) but with the integrated far field spectrum normalized by the integrated far field spectrum without 
the MNP. (f) as (c) but with the integrated far field spectrum normalized by the integrated far field spectrum in free space. 
The thin grey lines in (a),(b),(d),(e) indicate the LSP resonance. 



x-oriented and z-oriented dipoles, respectively, for both 
MNP sizes. Again the (invalid) dipole solutions are plot- 
ted for reference. The Lamb shifts at this height are quite 
large, giving a maximum frequency shift of |Aw| max /aj = 
7.9 x 1CP 3 for the 7-nm particle and |Aw| max /w = 1.28 x 
1CP 2 for the 20-nm particle in the x direction, and a max- 
imum frequency shift of |Aa>| max /u; = 2.02 x 10~ 2 for the 
7-nm particle and \A.u)\ max ./u) = 2.75x 10~ 2 for the 20-nm 
particle in the z direction. For comparison, at uj ~ 3 eV, 
an exciton linewidth of 7^ ~ 15mcV corresponds to 
7d/cj pa 5 x 10~ 3 , so the largest frequency shift in 
Fig- 121(f) is more than 5 times the exciton linewidth (even 
at room temperature), which, to our knowledge, is much 
larger than any previously reported result. For pho- 
tonic crystal systems 3 -!, \Aui\ max /uj «4x 10~ 5 has been 
reported, and for negative index metamaterial slaba 2 ^, 
|Aw| max /w «5x 10~ 4 has been predicted. 

B. Strong coupling regime and emitted spectrum 

Motivated by the significant enhancements seen in 
Fig. 03 we next study the nonpcrturbative strong- 
coupling regime, and calculate both the particle (or po- 
larization) spectrum and the spontaneous emission spec- 
trum of the field (Figs. [3] and H]). As mentioned above, 
for all calculations we use 7^ = 15 meV which corre- 
sponds to the decay of a typical QD exciton at room tem- 



perature 3 ^. Such a large decay would completely domi- 
nate semiconductor cavity systems, where the best (max- 
imum) vacuum Rabi splittings are around 0.1 meV— . For 
the 7 nm particle, with an x-oriented QD, the non-dipole 
result (i.e., not treating the MNP as a dipole) for P (lu) 
[Fig. 03 (a)] shows that there is a clear anticrossing, and 
the spectral location of strong coupling is evidently not 
located at the dipole LSP of 2.7915 eV; rather, it is much 
higher in energy at 2.9415 eV corresponding to the peak 
of the LDOS (~ 2.9489 eV). In contrast, the dipole result 
shows no indication of strong coupling. When looking at 
the far field spectrum [see Fig. 03(c)], it is more difficult 
to observe an anticrossing as the weighting provided by 
p nl causes the peaks to broaden and become more asym- 
metric. Additionally, there is a clear peak at 2.7885 eV 
which corresponds to the peak in p , this spectral peak 
is not observable in the particle spectra; this additional 
peak also shows up when using the dipole approximation 
but we emphasize that it is only due to photon propaga- 
tion from the MNP /emitter system to the detector and is 
extremely small. We highlight that these additional spec- 
tral features in the emission spectrum are quite different 
to a dielectric cavity system. 

From the calculations above, it is also clear that the 
predicted Lamb shifts are observable in the spectra as 
the exciton spectral peaks in both the particle and the 
far-field spectra are substantially shifted in energy. Fig- 
urcsO^b) and OHd) show the particle and light spectra for 
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the 20-nm MNP, and we observe many similar features 
to the 7-nm spectra; however the splitting between the 
peaks in the particle spectrum is notably larger for the 
20 nm particle compared to the 7 nm particle. In the far 
field spectra, we also see that there are significant quali- 
tative differences as the QD frequency is tuned; however 
the anticrossing region is observed at similar QD detun- 
ings. This shows that the strong coupling is clearly an 
observable effect, even with metal losses. 

Figure 5] shows similar results to those shown in Fig. [31 
but with a z-oricnted cxciton (dipolc). The larger LDOS 
for this polarization manifests in an increased Rabi split- 
ting in both particle and far field spectra, however the 
peaks at higher energies are more difficult to observe on 
this scale due to the smaller value of the non-local prop- 
agator in the higher frequency range [see Fig.^e)]. 

To further examine the anticrossing behavior of strong 
coupling, we have located the maxima in the particle 
spectra (red circles) and far field spectra (blue crosses) 
for various QD transition frequencies, and show these in 
Fig. [5] (a) for 7-nm particles, and Fig.[5](d) for 20-nm par- 
ticles using an x-oriented QD (similar results are found 
for a z-oriented dipole). In both the particle and light 
emission spectra, a clear anticrossing is observable indi- 
cating a vacuum Rabi splitting of around 2g m 79 meV 
for 7-nm particles and 2g rj 95 meV for 20-nm particles; 
note that 2g > 120 meV for both sizes of MNP for the z- 
oriented QD (not shown). Clearly this observation does 
not correspond to the location of the lowest-order dipolc 
mode (indicated by the thin grey lines in Fig. [3j but is 
due to the coupling to higher order plasmon modes. For 
the emitted light spectrum, the effects of propagation add 
additional peaks, however the vacuum Rabi splitting is 
well maintained even with non-radiative quenching. This 
finding is not at all clear unless one properly accounts for 
propagation to the detector. 

In Figs. OUa) and (d) we discern an additional peak 
in the far field at the location of the dipole mode which 
is due to light propagation (via p ni ), and in Fig. [5] (d) 
there is a fourth peak which is due to the dip located in 
p nl that only occurs for the 20-nm particle, but becomes 
too small to resolve after u>d ~ 2.85 eV. These spon- 
taneous emission spectra contain highly non-Lorentzian 
lincshapes as well as essential non-dipolar interaction ef- 
fects. Furthermore, any predictions of strong coupling 
with MNPs must include higher order mode coupling as 
they will dominate the dynamics before it is ever possi- 
ble to achieve strong coupling using the dipole mode (at 
least for our chosen parameters). 

Finally, we study the quenching effects in more detail. 
In Figs, in^b) and [5Je) we calculate the integration of 
the particle/far field spectrum as a function of u> d , and 
we normalize this to the integrated free-space value. We 
define the following integrated spectral quantities, Ip{uid) 
(integrated particle spectrum ) or Is{u<i) (integrated far 



field spectrum), which are computed as follows: 



Ip{ud) 



Is(ud) 



J °° P(u,u d ) doj 
Jo" S hom (uj,uJd) duj 



(17) 
(18) 



These integrals give the likelihood of detecting a pho- 
ton emitted by a QD, and the values in the vicinity 
of the MNP are normalized to the values that would 
be obtained from a QD in vacuum (for this particu- 
lar particle/detector geometry). We show MNP dipolc- 
approximation (dashed) and non-dipole (solid) results for 
a 7-nm MNP (orange-light line) and a 20-nm MNP (blue- 
dark line). For the integrated particle spectra, Ip, we 
see that, in terms of emitted flux, quenching is much 
more problematic for the non-dipolar result compared to 
the dipolar result. The region of greatest quenching is 
where the LDOS is peaked giving a maximum reduction 
to Ip w 0.3 in the region of the anti-crossing. Such an 
observation would lead one to believe that MNPs appear 
to absorb the majority of the emitted photons. However, 
in the far-field spectrum, Is, we see a dramatic increase 
in the relative number of photons detected for QDs lo- 
cated near the LSP of 30 (80) for 7 nm (20 nm). Even 
in the anticrossing region, the enhancement is « 2 — 3 
compared to a QD in free space. This enhancement in 
the integrated far field spectrum shows that even in the 
frequency region where photons appear to be dominated 
by non-radiative effects, the MNP compensates by acting 
as an antenna making the detection of far field radiation 
more efficient. 

To help further clarify the physics of metallic quench- 
ing, we also compare the MNP case with a metallic half 
space, where we calculate the Green function using a 
well known multilayer techniqu e 26 ! 38 . We initially ver- 
ify in Fig. [5jc) that a simple metallic half space suffers 
similar quenching to the MNP in the particle spectrum; 
however, it feels much more quenching in the far field 
spectrum [Fig. EJf)] with Is lower by about two orders 
of magnitude for both z-oriented (blue-dark line) and x- 
oriented (orange-light line) QDs. It is worth noting that 
the reduction of jd to values typical for QDs at cryogenic 
temperatures (~ [J-eV) results in much more dramatic 
quenching in both the particle spectrum and the far field 
spectrum for QDs coupled to MNPs; for example, using 
7d < 50 peY results in Is < 1 over the entire frequency 
range showing that the antenna effect of the MNP is un- 
able to overcome the quenching in the case of sharp QD 
lincwidths. 

For these quantum optical studies above, we have de- 
liberately chosen a rather large dipole moment (d = 24 
Debye) to enable the strong coupling regime. For smaller 
dipole sizes, e.g., with d = 12 Debye for the z-oriented 
dipole, or d = 8 Debye for a z-oriented dipole, we ob- 
tain qualitatively similar strong coupling results but with 
smaller vacuum Rabi splittings. There is also the poten- 
tial to see strong coupling with even lower QD dipolc 
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moments, if one uses MNPs with non-spherical shapes, 
e.g. cigar shapes^. For much lower dipole moments then 
the strong coupling effect of course vanishes, although 
dimer— ^ configurations may help to increase the LDOS 
to a sufficiently larger value. 

IV. DISCUSSION 

We briefly discuss some potential experimental config- 
uration for observing the effects presented above. There 
are several possible experimental scenarios that are likely 
within reach of current nano-fabrication techniques^. 
One example could involve spin coating colloidal QDs 
onto a substrate, locating the dots by correlating photo- 
luminescence data with atomic force microscopy (AFM) 
images and positioning the MNPs in the vicinity of the 
QD using the AFM tip as was done by Ratchford et al&; 
in their study, the relatively small QDs had an estimated 
dipole moment of around 5.3 Debye and a strong modifi- 
cation of the spontaneous emission rate was shown, along 
with a drastic reduction in blinking. A second possible 
method for probing the emission spectra of coupled QD- 
MNP systems could use an array of MNPs placed on a 
substrate and immersed in a solution of colloidal QDs. 
By illuminating with focussed (off-resonant) laser beams 
it is also possible to create efficient optical traps^— 
at which point the QD can be loaded into an excited 
state. Both of these proposals involve the use of sub- 
strates; of note our formalism enables the calculation of 
the far field spectrum in any inhomogeneous geometry as 
long as the Green function can be calculated; our initial 
results using the finite-difference time-domain technique 
[see Fig. Wijo)] for this simplified geometry can easily in- 
clude a substrate; or even more complicated geometries 
could be investigated such as MNPs coupled directly to 

I 



photonic crystal cavities containing QD o °' . In fact, 
recent experiments with QDs coupled with disordered 
metallic films on glass substrates are at a loss for the 
expected far field emission spectra^, further emphasiz- 
ing the usefulness of our technique. 

V. CONCLUSIONS 

We have presented a Green-function quantum optics 
approach to study quantum optical interactions between 
a dipole emitter and a single MNP. We began by exam- 
ining the properties of the classical Green function above 
a MNP within and beyond the dipole approximation and 
showed the dramatic effects of the higher order plasmon 
modes on the LDOS and photonic Lamb shifts. Going 
beyond the weak coupling approximation, we then ex- 
amined the particle spectrum and contrasted this with 
the far field (observable) light-emission spectrum of a 
QD coupled to the MNP. Using experimentally accessi- 
ble parameters, our non-perturbative light spectra show 
clear signatures of the strong coupling regime; the emit- 
ted spectrum was found to contain a triplet or quartet of 
resonances, highlighting the important role of light prop- 
agation to the detector. Finally, we also examined the 
role of quenching on the far field spectra, and compared 
the quenching to the case of a metal half space. Our tech- 
niques are quite general and can be extended to include 
an initial pump field, multiple MNPs and multiple QDs. 
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Appendix: Spherical Green function 

Given a sphere with permittivity e m and radius a, embedded in a homogeneous medium of permittivity e b , the 
scattered part of the Green function is given by 

G scatt (r,r')=-^££ E( 2 ~ 7 T , ■ (k b r)M% n (k b r') + R v N% n (k b r)N% n (k b r>)] , 

4tt n in + 1) (n + m)\ 

e,o n—1 m— v 

(A.l) 

where Rh/Rv are the centrifugal reflection coefficients corresponding to transverse electric/magnetic waves (TE/TM), 
and M™ n /N^ n are the vector functions corresponding to TE/TM waves and they have been separated into even and 
odd contributions. The values of R^ , Ry, are given by 

j-^h k m dT m T b k b dT b T m rX = ^ mTm ^ T b k b T b dT m ^ 

h III 

where 

n = i n (ha) , K t = h£> (kid) , (A. 3) 
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1 d(haj n (ha)) 
ba dha 



1 d(k l ah^ {ha)) 
kid dha 



(A.4) 



Here j n , hit are the spherical Bessel functions and spherical Hankel functions of the first kind respectively. The 
vector functions are defined as follows: 



M* mn (fcr) 



sin 6* 



d P m (cos / 

(fcr) P™ (cos 9) smm(j)0 - fcW (fcr) —2-4 ~ 



d0 



(A.5) 



777< 



-h^ (fcr)P™ (cosfl) cosm^ - (fcr) — 2_L i s in?^ 

) do 



(A.6) 



N^ n (fcr) = n ( n +l \ w (kr) P™ {cos 9) cos m# 
fcr 



f d(rh^(kr)) 



kr 



dr 



dP™ (cost 



cos mc 



m 
sin (9 



P™ (cos 9) sin mc; 



(A.7) 



N^„ (fcr) = n (rc+l) fe( i) (fcr) (cog ^ s . n 
fcr 



f d(rh^(kr)) 



(Zr 



dP™ (cosf 
~d9 



■ sin mc 



sin 9 



P™ (cos 9) cos m(f 



(A. 



Note that for our numerical calculation in this paper, a few simplifying assumptions can be made; we only consider the 
calculation to be along the z direction, x = x' = y = y' = 0, and we additionally assume that we are only calculating 
the LDOS (z = z'). This means that calculating in the 99 direction is equivalent to the cjxj) direction. This allows us 
to simplify Eqs. (fA~5]) - (|A~8]) to, 



M^(fcr) =hV (kr)^M^ 



(A.9) 



ML (fcr) =0, 



(A.10) 



N^„ (fcr) = n[n + l) h^ (kr) P™ (0) f , 



(A.ll) 
(A.12) 
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